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ABSTRACT 

In this study we investigate the formation and properties of prestellar and protosteUar cores using 
hydrodynamic, self-gravitating Adaptive Mesh Refinement simulations, comparing the cases where 
turbulence is continually driven and where it is allowed to decay. We model observations of these 
cores in the C^*0(2 — > 1), NH3(1, 1), and N2H+(1 — > 0) lines, and from the simulated observations 
we measure the linewidths of individual cores, the linewidths of the surrounding gas, and the motions 
of the cores relative to one another. Some of these distributions arc significantly different in the 
driven and decaying runs, making them potential diagnostics for determining whether the turbulence 
in observed star-forming clouds is driven or decaying. Comparing our simulations with observed cores 
in the Perseus and p Ophiuchus clouds shows reasonably good agreement between the observed and 
simulated core-to-core velocity dispersions for both the driven and decaying cases. However, we find 
that the linewidths through protosteUar cores in both simulations are too large compared to the 
observations. The disagreement is noticably worse for the decaying simulation, in which cores show 
highly supersonic infall signatures in their centers that decrease toward their edges, a pattern not seen 
in the observed regions. This result gives some support to the use of driven turbulence for modeling 
regions of star formation, but reaching a firm conclusion on the relative merits of driven or decaying 
turbulence will require more complete data on a larger sample of clouds as well as simulations that 
include magnetic fields, outflows, and thermal feedback from the protostars. 

Subject headings: ISM: clouds ~ kinematics and dynamics- stars:formation - methods: numerical — 
hydrodynamics - turbulence 



I. INTRODUCTION 

The origin of the stellar initial mass function (IMF) 
is one of the most important problems in astrophysics. 
Since the discovery of supersonic linewidths in star form- 
ing regions, understanding turbulence has been crucial 
for developing the theoretical framework for molecular 
cloud (MC) evolution, core formation, and the IMF. On- 
going debate in this field concerns whether the formation 
and destruction of MCs is dynamic and non-equilibrium 
(e.g. Elmegreen 2000; Hartmann 2001; Dib et al. 2007) 
or slow and quasi-equilibrium (Shu et al. 1987; Mc- 
Kee 1999; Krumholz et al. 2006b; Krumholz & Tan 
2007; Nakamura & Li 2007). The former mode would be 
characterized by transient turbulence, dissipating quickly 
on timescales comparable to the cloud lifetime so that 
GMCs are destroyed within ~ one dynamical time. The 
latter case corresponds to regenerated turbulence, per- 
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haps injected by the formation of the cloud, protostel- 
lar outflows, H II regions, external cloud shearing or su- 
pernova blastwaves, that is sufhciently strong to inhibit 
global gravitational collapse over many dynamical times. 
As shown by Offner et al. (2008) and Krumholz et al. 
(2005), the presence or absence of turbulent feedback di- 
rectly relates to the physical mechanism of star formation 
and determines whether stars form by the formation and 
collapse of discrete protosteUar cores (Padoan & Nord- 
lund 2002; McKee & Tan 2002) or competitive accretion 
(Bonnell at al. 2001). In the turbulent core model, the 
cloud remains near virial equilibrium on large scales and 
collapse occurs only locally in cores that are created and 
then mass-limited by the initial turbulent compressions. 
In the competitive accretion model, turbulence generates 
the initial overdensities, but without turbulent support, 
the cores are wandering accreting seeds, competing for 
gas from a reservoir, limited only by the size of the MC 
as a whole. 

There have been a number of recent observational pa- 
pers investigating starless and protosteUar core velocity 
dispersions, envelopes, and relative motions (Andre et al. 
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2007: Kirk ct al. 2007; Muench et al. 2007; Rosolowsky 
ct al. 2007; Walsh ct al. 2004), quantities that provide 
important clues about the core lifetimes and evolution, 
and about the turbulent state of the natal MC. All of 
these results, which include observations of a range of 
star forming regions in different tracers, indicate that 
observed low-mass cores have approximately sonic cen- 
tral velocity dispersions, at most transonic velocity dis- 
persions in their surrounding envelopes, and relative mo- 
tions that are slower than the virial velocity of the par- 
ent environment. Such results potentially contradict core 
properties mcasiired in simulations in collapsing clus- 
ters exhibiting competitive accretion (AylifFe et al. 2007: 
Klessen ct al. 2005; Tillcy & Pudritz 2004). 

In this paper we analyze the simulations described in 
Offner et al. (2008), which follow the evolution of an 
isothermal turbulent molecular cloud with and without 
continuous injection of energy to drive turbulent mo- 
tions. These simulations use the adaptive mesh refine- 
ment (AMR) code Orion (Truelove et al. 1998, Klein 
1999). The goal of our present work is to explore differ- 
ences between cores forming in these two environments 
and to provide predictions of their properties for observa- 
tional comparison. For this purpose, we simulate obser- 
vations of our cores using dust continuum and molec- 
ular lines, with realistic telescope resolutions. Unlike 
earlier comparisons of isothermal self-gravitating simu- 
lations with observations (Ayliffe et al. 2007; Klessen 
et al. 2005; Ballesteros-Paredes et al. 2003), we per- 
form more detailed radiative post-processing in order to 
simulate more accurately synthetic observations of our 
data. We also compare these observational measures 
for both driven and decaying turbulence, which has not 
previously been investigated. Keto & Field (2005) ob- 
tain post-processed simulated line profiles of several com- 
mon tracers modeled with a non-LTE radiative trans- 
fer code and find good agreement with observed iso- 
lated cores. However, their initial conditions arc sim- 
ple 1-D non-turbulent hydrostatic models and they halt 
the calculations when the central cores velocity exceeds 
the sound speed. Further, we report core-to-core cen- 
troid velocity dispersions of the simulated cores, which 
has not previously been studied in turbulent simulations. 
Work by Padoan et al. (2001) comparing observed large 
scale gas motions with 128'^ fixed-grid isothermal, non- 
self-gravitating, MHD simulations found good agreement 
with the gas centroid velocity dispersion-column density 
relation. In our higher resolution simulations, we instead 
focus on the smaller physical scales of self-gravitating 
cores and their observed properties, and we neglect the 
effects of MHD. 

In section 2, we describe our simulations in detail. Sec- 
tion 3 contains the methods of data analysis we use to 
simulate observations of our AMR data. In section 4, we 
present our results on the central core dispersions, rel- 
ative motions, and dispersions of the surrounding core 
envelopes. In section 5 we present quantitative compar- 
isons with observational data. Finally, section 6 contains 
our conclusions. 

2. SIMULATION PARAMETERS 

As described in Offner et al. (2008), our two simula- 
tions arc periodic boxes containing an isothermal, non- 
magnetized gas that is initially not self-gravitating. We 
first drive turbulent motions in the gas for two box cross- 



ing times, until the turbulence reaches statistical equi- 
librium, i.e. the power spectrum and probability density 
function shapes are constant in time. We adopt a 1-D 
Mach number of 4.9 (3-D Mach number of 8.5). At the 
time gravity is turned on, which we label t=0, our two 
simulations are identical. In one simulation energy injec- 
tion is halted and the turbulence gradually decays, while 
in the other turbulent driving is maintained so that the 
cloud remains in approximate virial equilibrium. The 
initial virial parameter is defined by 



GM 



= av 



~ 1.67, 



(1) 



where (Jid is the velocity dispersion, M is the cloud 
mass, and R = L/2 is the cloud radius. We use pe- 
riodic boundary conditions and 4 levels of refinement, 
which corresponds to an effective 2048^ base grid for an 
equal-resolution, fixed-grid calculation. 

Isothermal self-gravitating gas is scale free, so we give 
the key cloud properties as a function of fiducial values 
for the number density of hydrogen nuclei, nn, and gas 
temperature, T. It is then easy to scale the simulation 
results to the astrophysical region of interest. For the 
adopted values of the virial parameter and Mach number, 
the box length, mass, and 1-D velocity dispersion are 
given by 



L = 2.9 PC 

M = 865 Ti^/^fiH,V^ Mq, 
c7iD = 0.9 Ti^/^ km s"^ , 



tff = 1.37 njj^V^ Myr, 



(2) 
(3) 
(4) 
(5) 



where we have also listed the free-fall time for the gas in 
the box for completeness. 

These equations are normalized to a temperature Ti = 
T/10 K and average hydrogen nuclei number density 
'^H.s = "■h/(1 X lO'' cm^''). For the remainder of this pa- 
per, all results will be given assuming the fiducial scaling 
values of riH = 1-1 x 10^ cm~^ and T = 10 K (Perseus) 
or fiH = 2.0 X 10^ cm"^ and T = 20 K (p Ophiuchus; see 
§5) and assuming a mean particle mass of /i=2.33mH- 
These conditions place p Ophiuchus slightly above the 
observed lincwidth-size relation (Solomon et al. 1987; 
Heyer & Brunt 2004): 



(TiD = 0.5 



L 



l.Opc 



0.5 



km s 



(6) 



where L is the cloud length (we assume that Perseus lies 

on this relation-see §5 below). 

Note that this relation differs somewhat from the re- 
lation given by Heyer & Brunt (2004) since the length 
scale determined from a Principal Component Analysis 
is smaller than the actual size of the region being sam- 
pled (see McKee & Ostriker 2007). These parameters 
may be adjusted to different conditions using equations 
(2)-(5). However, once we simulate an observation of the 
data for a given tracer, the scaling is fixed. Using these 
units, the minimum cell size is ^ 90 AU and 280 AU for 
p Ophiuchus and Perseus, respectively. 

In the simulations, we introduce sink particles in col- 
lapsing regions that violate the Jeans condition (Tru- 
elove et al. 1997) at the finest AMR level (Krumholz et 
al. 2004), where we adopt a Jeans number of J = 0.25. 
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Cores that contain sink particles are analogous to ob- 
served protostellar cores, which contain a central source, 
while cores without sink particles can be considered 
prestellar. This distinction is an important one in some 
cases and we discuss some differences in the two simula- 
tions in §4. Note that due to our resolution and neglect of 
protostellar outflows, the sink particles represent a mass 
upper limit for any potentially forming protostar. 

3. ANALYSIS 

Since our goal in this paper is to contrast the simu- 
lations and compare them with observations, we must 
attempt to replicate an observer's view of our simula- 
tion. Observations of core kinematics, such as those of 
Andre et al. (2007, henceforth A07), Kirk et al. (2007, 
henceforth K07), and Rosolowsky et al. (2007, hence- 
forth RG7), generally trace the gas mass using dust con- 
tinuum data and obtain velocity information by observ- 
ing the same region in one or more molecular tracers. 
We process our simulations using a rough approximation 
of these techniques as follows. First, we select a fiducial 
cloud distance of either 125 pc, corresponding to the dis- 
tance to the Ophiuchus star-forming MC, or 260 pc for 
comparisons with the Perseus MC. Second, we select an 
appropriate telescope resolution of 26" or 31" FWHM, 
corresponding to 0.02 pc and 0.04 pc at our adopted dis- 
tances, and approximate the telescope beam as Gaussian 
in shape. We perform all line fits assuming 0.047 km 
s^^ velocity resolution per channel. Increasing the ve- 
locity resolution further has little effect on the line fits. 
For simplicity we adopt the same resolution for obser- 
vations in dust continuum and in all molecular tracers. 
Our fiducial resolution is typical of observations of core 
kinematics (e.g. A07, K07,'r07). 

For the dust continuum observations, since our gas and 
dust are isothermal and the simulation domain is every- 
where optically thin at typical observing wavelengths of 
~ 1 mm, the dust intensity emerging from a given pixel is 
simply proportional to the column density in that pixel. 
We therefore define a dust continuum map by comput- 
ing the column density and convolving the resulting map 
with the beam. To avoid introducing unnecessary and 
artificial complications, we neglect observational uncer- 
tainties in the conversion from an observed dust contin- 
uum intensity to a column density, and assume that the 
column density can be reconstructed accurately except 
for beam smearing effects. We identify cores by finding 
the local maxima directly from the column density data. 
In the analysis, we consider only local maxima with peak 
columns greater than 0.1 of the global maximum column 
of the smeared data. This cutoff corresponds to ^ twice 
the mean smeared column density. 

To model molecular line observations, we choose three 
representative lines, the J = 2 ^ 1 transition of C^^O, 
J = 1 ^ transition of N2H+, and the NH3(1, 1) tran- 
sition, which have critical densities of 4.7 x 10^ cm~^, 
6.2 X 10^ cm^'', and 1.9 x lO'' cm^"', respectively. (For 
this calculation and all the others presented in this paper, 
we use molecular data taken from the Leiden Atomic and 
Molecular Database^, Schoier et al. 2005). These fines 
are often used in observing core kinematics because they 
span a range of densities and, with the exception of C^^O 
along the densest sightlines, are generally optically thin 

^ See http://www.strw.leidenuniv.nl/~moldata 



in low mass star-forming regions. We discuss the issue of 
optical depths in more detail in § 5.2. 

We generate a position-position-velocity (PPV) data 
cube from our simulations in each of these lines us- 
ing a two step procedure, which combines elements of 
Krumholz et al. (2007a) and Krumholz et al. (2007b). 
The first step is to compute the emissivity as a function of 
density. Since, as we shall see, the density-dependence 
of the molecular emission has important consequences, 
we cannot assume that these species are in local ther- 
modynamic equilibrium (LTE). Instead, we assume that 
the gas is in statistical equilibrium, that it is optically 
thin, and that radiative pumping by lino photons is neg- 
ligible. Note that the advection time of the gas is large 
compared to the molecular coUisional and radiative time 
scales, which are on the order of a few years for the mean 
density of our simulations. Thus, the gas reaches statis- 
tical equilibrium essentially instantaneously relative to 
the gas motion. CoUisional excitation dominates over ra- 
diative excitation or de-excitation by line photons along 
lines of sight where the density is above the transition 
critical density. Since we are particularly interested in 
the high density regions of the cores, we need not con- 
sider radiative pumping in our analysis. However, we do 
include radiative excitation and de-excitation due to the 
cosmic microwave background, since this can be signifi- 
cant for lines at very low frequencies such as NH3(1, 1). 

For a molecule like C^^O with no hyperfine structure, 
under these approximations the fraction fi of molecules 
of a given species in bound state i is given by the equa- 
tions of statistical equilibrium 

y^(^H2gji + ^ji + BjiIcMB)fj 



^('^Haftfc + ^ik + BiklcMn) 



fi 



(7) 



where rin.^ is the molecular hydrogen number density, qij 
is the collision rate for transitions from state i to state 
j, A and B are the Einstein coefficients for this transi- 
tion, and /cMB is the intensity of the cosmic microwave 
background radiation field (which is simply the Planck 
function for a 2.73 K blackbody) evaluated at the tran- 
sition frequency. In this expression we adopt the con- 
vention that the summations run over all bound states, 
the spontaneous emission coefficient Aij = for i < j, 
that Bij is the stimulated emission coefficient for i > j, 
the absorption coefficient for i < j, and is zero for i = j, 
and that = for i = j. For molecules with hyperfine 
structure, we show in Appendix A that with some ad- 
ditional approximations equation (7) continues to hold 
provided that the rate coefficients Qij ^ Aiji and B-ij are 
understood as statistically- weighted sums over all the hy- 
perfine sublevels of states i and j. 

For molecules without hyperfine structure, the net 
emission minus absorption of the background CMB pro- 
duced by a parcel of gas along the line of sight is then 
given by 



3ij - Xijic 



MB 



hVij 



X [fi{-^ii + BijlcMs) — fjBjilcMB], 



(9) 



4 



where Xij is the extinction of the CMB due to resonant 
absorption, I'ij is the transition frequency, X is the abun- 
dance of the species in question relative to hydrogen nu- 
clei, and nn is the number density of hydrogen nuclei. 
Physically, this quantity represents the net radiation in- 
tensity added by transitioning molecules over and above 
what one would see at that frequency due to the CMB 
alone, under the assumption that the line is sufficiently 
thin that the CMB dominates the intensity at that fre- 
quency. It is the intensity one will observe in a line after 
subtracting off the continuum. In the case of a molecule 
with hyperfine structure, under the standard assTimption 
that the hyperfine sublevels are populated in proportion 
to their statistical weight (sec Appendix A), the inten- 
sity produced by a single transition from level i, hyperfine 
sublevel a to level j, hyperfine sublevel /3 is given by 

3iajf3 — Xiajp^CMB ~ 



-Xnii 



fi^^{Aiajl3 + Biajj3lcMB) — fjBjj3iaIcMB 

9i 9j 



(10) 



where gia is the statistical weight of hyperfine sublevel 
Q^j 9i = X^a 9i is the summed statistical weight of all 
the hyperfine sublevels of state i, and the combination 
of subscripts iaj/3 indicates the frequency or radiative 
coefficient for transitions from level i, hyperfine sublevel 
a to level j, hyperfine sublevel (3. If one neglects the 
very small differences in frequency between the different 



hyperfine transitions (i.e. 



takes 



indepen- 



dent of a and /?) and sums equation (10) over hyper- 
fine substates a and /?, then it immediately reduces to 
equation (9) provided that the rate coefficients are un- 
derstood to be statistically-weighted sums of the individ- 
ual hyperfine transition coefficients (per equations A6 - 
A8). Thus equation (9) gives the total intensity summed 
over all hyperfine components. In either the presence or 
absence of hyperfine splitting, to compute the intensity 
from our simulations, we solve the system of equations 
(7)-(8) for our fiducial temperature T for a wide range 
of molecular densities tihs and tabulate the quantities 
Uij - XijIcMB)/X or {jiaj/} - Xiaj/jlcMB)/ X as a func- 
tion of riHa • 

The second step to generate the PPV cube from the 
simulation data is to compute the emergent intensity in 
each pixel in each velocity channel using our tabulated 
net emission function. The specific emissivity minus ab- 
sorption of the gas at a frequency v is (jy —XijIcMB)4>{^) 
or {jiajp — Xiaji3lcMB)4>{v), in the absence or presence 
of hyperfine splitting, where (j){v) is the line shape func- 
tion. To determine 4'{v)^ we assume that the molecules 
in each cell have a Maxwellian velocity distribution with 
dispersion (j„ = y^fcsT/m, where m is the mass of the 
emitting molecule. For this velocity distribution, the line 
shape function for a fluid with bulk velocity v is 



(l>{vohs\v) 



1 



exp 



2al 



(11) 



where an observation at velocity Vohs is understood to 

mean an observation at frequency v = [1 — Vobs/cji'ij 
and where ai, = [a^/cjvij. For optically thin emission 
with no hyperfine structure at an observed velocity Wobs, 
a cell of length Aa; contributes a specific intensity above 
the continuum of 



where jij and Xij are functions of the cell density nn and 
<?i'(^'obs; v) is a function of the cell velocity v. The intensity 
averaged over a velocity channel that covers velocities in 
the range vq < Vohs < v\ is 



chan 



Hi] ^ XijIcMB) 



cAx 



erff^' 



(13) 



= Uij - XijIcMB)Ax(j){Vohs;v) 



(12) 



We compute the channel-averaged specific intensity along 
each line of sight by summing (/i/)^han '^^^^ cells, 
each with its own velocity v, along the line of sight. The 
final step in constructing our PPV data cube is that we 
take the summed intensity computed in this way and 
smear each velocity channel using our Gaussian beam. 

In the case of molecules with hyperfine structure, the 
equations are identical except that the subscripts ij are 
replaced by iajf3, and we note that, since the hyper- 
fine components are closely spaced in frequency, multiple 
components may contribute significant intensity at the 
same frequency. However, in the observations to which 
we wish to compare our simulations, kinematic informa- 
tion is generally obtained by fitting one or more well- 
separated individual hyperfine components (e.g. A07, 
K07, although see ROT, who use a more complex pro- 
cedure). Thus, in practice it is generally not necessary 
for our purposes to consider more than a single hyperfine 
component. For optically thin emission in hyperfine com- 
ponents with no significant line overlap, this means that 
the procedures for molecules with and without hyperfine 
splitting are the same. 

Our procedure determines the emission only up to the 
unknown abundance X, which in reality will depend on 
the emitting species and on the density and temperature, 
and probably also the thermal and density history, of a 
given fluid element. For example, observations show that 
in the densest cold regions CO and its isotopomers will 
be depleted, while the abundance of N2H+ stays roughly 
constant (Tafalla et al. 2004a, b). In order to approxi- 
mate this effect, we adopt a simple depletion model for 
each of the chemical species that we simulate. For C^^O, 
we assume an abundance of AT = 10~^ molecules per H2 
molecule with depletion occurring at riH2 = 5 x 10^ cm~^ 
(Tafalla et al. 2004a). For N2H+, we adopt X = 10-^° 
with depletion at riHs = 5 x lO'' cm~^ (K07; Tafalla et al. 
2002). Although depletion in nitrogenous species is not 
generally observed, it is assumed that N2 begins to dis- 
appear at number densities riHs > 10^ cm~^ (Walmsley 
et al. 2003). For the NH3 measurements we compare to 
in Perseus, we use X = 10~^ (Rosolowsky, private com- 
munication) with assumed depletion at the same density 
as N2H+. 

We use these procedures to produce dust continuum / 
column density maps and PPV cubes for each of our three 
molecular lines. To increase our statistics, we generate 
data sets for each cardinal direction a.t t = tg, and we 
treat the three orientations as independent observations. 
Figure 1 shows a dust continuum map in one particular 
orientation. 

4. RESULTS 

In the decaying simulation, at Ifff we identify a total 
of 109 cores, 54 of which can be considered protostellar 
due to the presence of a sink particle within 0.1 pc of 
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the core center. In the driven simulation, we find 214 
cores, 92 of which are protostellar. A large central point 
mass can have a significant effect on the core gas mo- 
tion, so we separate out the 'starless' cores for compari- 
son. The relative number of starless cores to protostellar 
cores varies from star-forming cloud to cloud depending 
upon the advancement of star-formation in the region. 
The ratios of prestellar to protostellar cores that we find 
in our simulations are similar to the ratios observed in 
Perseus and Ophiuchus (Young et al. 2006; Enoch et 
al. 2006). In these simulations, the larger number of 
cores in the driven run is significant because the ongoing 
turbulence creates more new condensations, which also 
collapse more slowly. 

For the sake of clarity, we will refer to the centroid 
velocities of the cores as the "first moments" and the ve- 
locity dispersions through the core centers as the "second 
moments." Thus in the following sections we will describe 
the measured distributions of the first and second mo- 
ments and report the dispersion of the first moments (i.e. 
the core-to-core velocity dispersion). We define transonic 
velocities as those falling in the range Cs < ct < 2cs, while 
supersonic dispersions have u > 2cs. 

4.1. Central Velocity Dispersions 

In this section, we investigate the distribution of sec- 
ond moments (central non-thermal velocity dispersions 
through the core centers) in N2H+, a measure that is 
uscfiil for determining the level of turbulence and infall 
motion in the core. The total dispersion along the line 
of sight is given by 



where ut = \Jk'QT jra and ctatt is the non-thermal com- 
ponent that we discuss here. 

We compute ctlos in the simulations by fitting a Gaus- 
sian to the spectrum through the core center and then 
deriving the second moment, ctnt, from equation 14. Ta- 
ble 1 gives the median and means of unt /cs, and we plot 
the total distribution in figure 2 and the prestellar and 
protostellar distributions in figures 3 and 4, respectively. 
The core populations appear fairly similar in the two 
simulations, although there is evidence of the increased 
turbulence in the driven simulation. Since the cores are 
created by turbulent compressions in both environments, 
at early times they should have similar second moments. 
However, at late times, as the cores collapse and form 
protostars the distributions are more dissimilar. Indeed, 
from figure 4 we can see that the protostellar distribu- 
tions are much broader and less peaked than the prestel- 
lar ones. The decaying protostellar core population has 
almost twice as many cores in the tail (ctnt > 4cs) of the 
distribution, while the protostellar driven population is 
dominated by cores with ctnt < 4cs. 

To better characterize the differences between the two 
simulations, we perform a Kolmogorov-Smirnov (KS) 
test comparing each of the core distributions. The KS 
statistic gives 1 minus the confidence level at which the 
null hypothesis that the two samples were drawn from 
the same underlying distribution can be ruled out, e.g. 
a KS statistic of 0.01 means that we can reject the hy- 
pothesis that the two samples were drawn from the same 
distribution at the 99% confidence level. We find that the 
net driven and decaying velocity dispersion populations 



have a KS statistic of 18%, meaning that we can rule out 
the hypothesis that they were drawn from the same pop- 
ulation only with 82% confidence. Individually, there is 
large disagreement in both the protostellar populations 
(4 X 10~^%) and prestellar core populations (2%). 

The difference between the protostellar populations in 
the two simulations is associated with the mass differ- 
ences between the sink particles: The decaying simula- 
tion has a median sink mass that is approximately twice 
that of the driven simulation and correspondingly larger 
accretion rates that are associated with higher velocity 
dispersions. 

4.2. Gore Envelopes and Surroundings 

The velocity dispersions of the gas surrounding the cen- 
tral column density maxima yield information about the 

relative motion between core and envelope, and may also 
reveal the presence of shocks or strong infall that could 
limit core boundaries. Typically, observers find only 
small differences in velocity between the core and the 
surrounding gas envelope, which rules out dynamical pic- 
tures of core accretion in which protostars may strongly 
gravitationally interact with their neighbors (K07). In 
addition, although shocks arc postulated to be the origin 
of the original density compression, close observations 
have not revealed strong confining shocks surrounding 
the cores. Generally, our simulations produce prestel- 
lar cores that agree with the expectations from observa- 
tions. However, the decaying protostellar cores exhibit 
supersonic internal velocities that are not observed in the 
star-forming regions we compare with. 

In order to compare the two environments observed 
with three common tracers, C^^O, N2H+, and NH3, we 
calculate the velocity dispersion through each pixel along 
the line of sight. Figures 5 and 6 show the velocity dis- 
persion of each pixel in the vicinity of a single prestellar 
and protostellar core for decaying turbulence, which rep- 
resent typical examples of each type from our sample, 
overlaid with contours of integrated intensity. The large 
number of cores in our sample makes comparing the pop- 
ulations by eye on an individual basis difficult. In order 
to consolidate the data sets for each environment, we bin 
the pixels by radial distance from the core center. We de- 
fine 20 logarithmic bins that range from 0.005 to 0.1 pc 
in projected distance from the core center and then av- 
erage together the velocity dispersions of all pixels that 
fall into a given bin, including all prestellar or protostel- 
lar cores in each case. The result is a single 'averaged' 
core for each tracer and environment. We have plotted 
this averaged velocity dispersion as a function of distance 
from core center in figures 7 and 8 for starless and proto- 
stellar cores, respectively. There are several interesting 
points that may be noted from these plots. 

First, gas sampled by low density tracers (e.g. C^^O) 
around prestellar cores has a higher velocity dispersion 
than that sampled by higher density tracers. This is rea- 
sonable given that the lower-density gas is further from 
the core center and generally more turbulent. Before 
collapse ensiles, the cores have typically not developed 
strong high density peaks as is evident in figure 5. This 
difference between lower and higher density tracers has 
been frequently exploited observationally to distinguish 
between the dense core and surrounding envelope (e.g. 
K07; Walsh et al. 2004). 

Second, figure 7 shows that the starless cores forming 
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Fig. 1. — The images show the decaying (left) and driven (right) log column densities (g cm~^) 'observed' at a distance of 260 pc with 
beam size of 31" . 



TABLE 1 

Central velocity dispersion median and mean for the two environments and core types at l.Oifj in N2H+ normalized to the 

CONDITIONS IN PeRSEUS. 



Decaying Driven 
All Prestellar Protostellar All Prestellar Protostellar 



A^coros 109 55 54 
Median (tnt/cs 1.0 0.6 2.9 
Mean ct^t /cg 2^2 06 3^ 

in the driven simulation tend to have a higher average 
velocity dispersion than those in the decaying simulation. 
This is mainly apparent in the tracer C^^O, which traces 
the more turbulent core envelope. 

Most importantly, the average prestellar velocity dis- 
persion for both cases and for all tracers are approxi- 
mately sonic. Even the lowest density tracer, C^^O, re- 
mains, on average, below 2cs for the range of column 
densities in the core neighborhood. 

Finally, we note that there is only a small increase in 
the dispersion with increasing radius. This is consistent 
with observations by Barranco & Goodman (1998) and 
Goodman et al. (1998) who find that the velocity disper- 
sion of the cores on the scale of '--^ 0.1 pc is approximately 
constant, with a small increase near the edge of this re- 
gion of "coherence." The magnitude of the dispersion 
suggests that the starless cores forming in a turbulent 
medium are not strongly confined by shocks in the range 
of densities that are traced by observers. 

In contrast, some of these conclusions do not hold for 
protostellar cores, when strong infall occurs. As shown 
in figure 8, protostellar cores exhibit significantly higher 
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1.1 0.9 2.1 

1.8 1.2 2.7 



average velocity dispersions than the prestellar counter- 
parts. The tracers of the protostellar cores behave dif- 
ferently as well. Due to the strong infall, which occurs 
in the densest gas, the higher density tracers, N2H+ and 
NH3, show higher velocity dispersions than the C^^O, 
which indicates that the lower density envelope remains 
transonic. 

There is also clearly a significant difference between the 
protostellar cores in the two environments. Those cores 
in the driven environment have transonic to slightly su- 
personic velocity dispersions in all tracers that do not 
vary significantly with distance from the core center, 
which is consistent with the coherent core structure ob- 
served. This indicates that the cores still have residual 
turbulent pressure support at a global freefall time and 
collapse more slowly. However, the protostellar cores in 
the decaying turbulence environment, lacking this sup- 
port, have shorter lifetimes and proceed more quickly 
to collapse and develop much higher, supersonic, cen- 
tral velocity dispersions in N2H+ and NH3 as the cloud 
gas infalls to the high density regions. At large radii 
however, the velocity dispersion of the protostellar cores 
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Fig. 2. — Fraction / of all cores binned as a function of second moments (non-thermal velocity dispersion), (Tnti for a simulated observation 
of Perseus using N2H+. The distribution on the left shows the cores in the decaying turbulence enviornment, while the distribution on the 
right gives the cores in the driven turbulence enviornment. 
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Fk;. 3. — Fraction / of starless cores binned as a function of second moments (non-thermal velocity dispersion), (Tnt, for a simulated 
observation of Perseus using N2H^. The distribution on the left siiows those cores in the decaying turbulence enviornment, while the 
distribution on the right gives the cores in the driven turbulence enviornment. 



in the decaying enviroment matches the velocity disper- 
sion of cores in the driven environment. A similar time- 
dependent trend is obtained in decaying simulations by 
Ayliffe et al. (2007). 

In summary, prestellar cores forming in driven turbu- 
lence have average dispersions of ^ 1.5cs in all tracers, 
and this dispersion is either flat or slowly decreasing with 
increasing radius. In contrast, cores in decaying tur- 
bulence show small (ctnt < l-Ocg), flat dispersions for 
prestellar cores, but large and radially decreasing disper- 



sions for protostcUar cores. This is most likely due to 
infall of unbound gas from large distances at late times, 
which is a signature of competitive accretion. We do 
not observe this in the driven run because the cloud gas 
dispersion is too high for Bondi-Hoyle accretion to be 
efficient over large distances (Krumholz et al. 2006a). 

The dispersions we obtain for the cores and their sur- 
rounding envelopes are somewhat dissimilar to those ob- 
tained by Klessen et al. (2005) in SPH simulations. As 
we do, Klessen et al. investigate the velocity dispersions 
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Fig. 4. — Fraction / of protostellax cores binned as a function of second moments (non-thermal velocity dispersion), ctnti for a simulated 
observation of Perseus using N2H+. The distribution on the left shows the cores in the decaying turbulence enviornment, while the 
distribution on the right gives the cores in the driven turbulence enviornment. 



of cores forming in an isothermal, large scale driven tur- 
bulent environment. In their study, they derive clump 
properties when only 5% of the mass is in cores or at 
^ OAtjf, a much earlier time than wc use. However, 
even for prestellar cores with driving, they frequently 
find strong supersonic shocks with ctlos ~ 3 — 5cs bound- 
ing the cores, which is thus far not supported by obser- 
vations. In lieu of a simulated observation, they use a 
column density cutoff to make the dispersion estimates. 
We find that wc obtain higher velocity dispersions cal- 
culating the velocity dispersion directly as Klessen et al. 
do rather than fitting the line profile in the manner of 
observers. The reason for the difl[erence is that in some 
cases the spectra resemble a fairly narrow peak, which 
is well fit by a Gaussian, surrounded by a much broader 
base around the 10% level. The magnitude of this extra 
spread is reduced substantially at the higher densities as 
traced by N2H+, and it is likely neglected in the fits per- 
formed by observers due to the inherent low-level noise in 
the actual spectra. Another possibility for the difference 
is the difficulties of SPH in rendering shocks and instabil- 
ities, in particular shear flow instabilities (Agertz et al. 
2007) that are likely to be present in any compressible 
turbulent simulation and may seriously affect accuracy. 
However, the extent that this may contribute to the high 
dispersions found by Klessen et al. is unclear. 

4.3. Relative Motions 

Observers frequently evaluate an intensity-weighted 
mean velocity, or first moment, along the line of sight 

through the core center. While the second moments are 
indicative of infall motions, the first moments represent 
the net core advection. The dispersion of the first mo- 
ments indicates how much the cores move relative to one 
another. Observations find that the dispersion of first 
moments is generally smaller than the velocity disper- 
sion of gas that is not in cores, although how much so 
varies from region to region. For example, A07 conclude 



that the first moment dispersion is sub-virial by a fac- 
tor of ^ 4 in p Ophiuchus. K07 find that first moment 
dispersion of starless cores in Perseus is sub-virial by a 
factor of ^ 2, which does not rule out virialization. 

In order to get an unbiased distribution for comparison, 
it is necessary to subtract out any large gradients in the 
sample of first moments. Thus, for each region we fit 
V = Vo + VV ' X clS cl function of position, x. Generally, 
this turns out to be a fairly small correction, but the 
net effect is to decrease the dispersion of first moments 
relative to the gas 

We plot the distribution of first moments for all cores 
in both environments in figure 9, and we plot the dis- 
tributions for prestellar and protostellar cores separately 
in figure 10. In these, we normalize to the "measured" 
gas dispersion and correct for the velocity gradient in the 
box. The dashed line is a Gaussian with the same dis- 
persion as the core distribution. For reference, we also 
plot a Gaussian with the gas dispersion. Note that in the 
driven simulation the dispersion inferred from virial ar- 
guments and the time-dependent gas dispersion are the 
same, because by definition we fix the total kinetic en- 
ergy to maintain virial balance. However, for the de- 
caying simulation, the time-dependent gas dispersion is 
lower than would be derived from a virial argument using 
the total gas mass and cloud size. 

Again, we use KS tests to characterize similarity in 
the populations, which we report in Table 2. A KS test 
indicates that driven and decaying distributions of the 
net first moments agree with 56% confidence, while the 
prestellar and protostellar core first moments agree with 
40% and 13% confidence. This is significant enough to 
imply that the early core motions arc not widely differ- 
ent in the two environments, with the largest difference 
occurring between the protostellar first moments. Com- 
paring these distributions with a Gaussian dispersion at 
the gas dispersion yields good agreement for the distribu- 
tions of the prestellar driven cores (54% confidence) and 
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Fig. 5. — The upper plot gives average velocity dispersion as a function of radius for a single decaying starless core at Itg. The images 
below show a simulated observation in C^*0 (left) and N2H"'" (right). Contours indicate integrated intensity where each contour is a 10% 
linear change from the peak specific intensity in that tracer. The color scale shows velocity dispersion, (tnt/cs and the circle indicates the 
FWHM beam size. 



protostellar decaying cores (56%), but low agreement for 
the other distributions. In general, low agreement may 
be because the first moment distributions, although hav- 
ing a similar dispersion to the gas in some cases, are not 
well represented by a Gaussian distribution. 

In Table 3, we list the first-moment dispersions, both 
corrected and uncorrected for large linear gradients. We 
find that the corrected net core dispersion for the driven 
and decaying cores are both sub-virial relative to the gas 
dispersion. Previous simulations have shown that the 
dispersion of first moments becomes sub-virial towards 
higher gas densities (Padoan et al. 2000), so the result is 
not unexpected. One interesting difference between the 
simulations is that the decaying protostellar cores are ap- 
proximately virial, while the prestellar driven cores are 
approximately virial. The former suggests that as the 
cloud loses turbulent support and tends toward global 



collapse, that either the core interactions increase or that 
the cores retain some memory of their natal gas disper- 
sion. The inertia of the cores implies that their velocity 
dispersions will tend to decay more slowly than that of 
the whole. This is a potentially testable signa- 

ture of the competitive accretion model (Bonnell et al. 
2001). In the latter case, the prestellar cores may still 
be forming out of the shocking gas and hence may still 
have similar motions. In general, the sub-virial disper- 
sion of the cores may imply that they are not scattering 
sufficiently frequently to virialize within the formation 
timescale. Elmegreen (2007) reasons that if cores form at 
the intersection of two colliding shocks, then their initial 
dispersion should be on average less than the gas disper- 
sion. Overall, our results imply that the forming cores 
are at least somewhat sensitive to the actual dispersion 
of the natal gas. 
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Fig. 6. — The upper plot gives average velocity dispersion as a function of radius for a single decaying protostellar core at Itg. The 
images below show a simulated observation in O (left) and N2H+ (right). Contours indicate integrated intensity where each contour 
is a 10% linear change from the peak specific intensity in that tracer. The color scale shows velocity dispersion, (tnt/cs and the circle 
indicates the FWHM beam size. 



5. OBSERVATIONAL COMPARISONS 

5.1. Scaling to Observed Regions 

In this section, we compare our simulated observations 
with three selections of cores observed in three standard 
molecular tracers in two different low-mass star-forming 
regions, p Ophiuchus (primarily L1688) and the Perseus 
Molecular Cloud. This comparison cannot be precise for 
several reasons: First, the cloud is isolated, whereas our 
simulation is a periodic box; second, we are using a sin- 
gle simulation with given values of the virial parameter 
and the Mach number to compare with clouds that have 
somewhat different values of each of these parameters; 
and, finally, our simulation is isothermal, whereas the 
temperature is observed to vary in the clouds. Further- 
more, the actual cloud is magnetized, whereas our sim- 
ulation is purely hydrodynamic. A variety of possible 



comparison strategies is possible. We have chosen to use 
the same mean density in the box as in the cloud, and 
to make the simulation temperature agree approximately 
with the typical temperature observed in the cloud cores. 
The size and mass of the simulation box then follow from 
equations (2) and (3). With this approach, the Jeans 
mass will be about the same in the simulation and in 
the cloud, but the size and mass of the overall cloud will 
generally differ between the two. 

A07 observed 41 starless cores in p Ophiuchus and 
made maps of 26 of them using the tracer N2H+ (J = 
1^0), which are clustered in a region of area 1.1 pc^. 
The total gas mass in this region with extinction greater 
than 15 magnitudes is estimated to be ^ 615 Mq (Enoch, 
private communication; Enoch et al. 2007) with peak col- 
umn densities of 7Vh2=1-8x 10^^ cm^^ (Motte & Andre 
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Fig. 7. — The figures show the averaged dispersion of the prestellar cores binned over distance from the central core, where D denotes 
driven and U denotes undriven turbulence. 
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Flc;. 8. — The figures show the averaged dispersion of only the protostellar cores binned over distance from the central core, where D 

denotes driven and U denotes undriven turbulence. 



1998). The star- forming area of p Ophiuchus is roughly 
circular with radius R ^ 0.6 pc; the mean density and 
column density are therefore nn — 2 x 10"' cm""^ and 
A^H = 5 X 10 cm~^. As discussed above, we adopt 
this density for our simulation. To fix the temperature, 
we first note that dust temperatures in the pre-stellar 
cores range from 12-20 K (A07). On the scale of the en- 
tire L1688 cloud, the temperatures as measured by ^^CO 
and ^^CO lines are 29 K and 21 K, respectively (Loren 
1989a; in his notation, this region is R22). We therefore 
adopt T = 20 K for the simulation. Equations (2) and 
(3) give L = 0.9 pc and M = 550 Mq for the simula- 



tion box, comparable to, although somewhat less than, 
the ol)servcd values. The total velocity dispersion mea- 
sured from the '^'^CO line is 1.06 km s~^ (Loren 1989b), 
which lies above the standard linewidth-size relation (eq. 
6). The corresponding ID Mach number is A4id = 3.9, 
slightly less than the value 4.9 in the simulation. The 
virial parameter of the cloud is 1.25, also slightly less 
than the simulation value of 1.67. 

For the Perseus MC, K07 report central velocity dis- 
persions and ccntroid velocities measured from C^^O 
and N2H+ pointings for 59 prestellar and 41 protostellar 
cores. R07, also making pointed observations of Perseus, 
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Fig. 9. — Fraction / of all cores binned as a function of first moments, Vcent, for a simulated observation using N2H+ normalized to 
the large-scale gas dispersion. Vg a,t t = tg. The distribution on the left shows the cores in the decaying turbulence environment, while 
the distribution on the right gives the cores in the driven turbulence environment. The dashed line is a Gaussian with the same dispersion 
as the data while the dot-dashed line is a Gaussian with the gas velocity dispersion (Vg = 2.2cs, = 4.9cs, for the decaying and driven 
simulations, respectively). 




Fk;. 10. — Fraction / of prcstcUar cores (top) and protostcUar cores (bottom) binned as a function of first moments, Vcent > for a simulated 
observation using N2H+ normalized to the large-scale gas dispersion, Vg. The distribution on the left shows those cores in the decaying 
turbulence enviornment, while the distribution on the right gives the cores in the driven turbulence enviornment. 
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TABLE 2 

KS STATISTICS FOR THE DRIVEN AND DECAYING CORE FIRST 
MOMENTS (CENTROID VELOCITIES) CORRECTED FOR LARGE 
VELOCITY GRADIENTS AND THE GAS. 



from state i to state j is 







D: All 


D: Starless 


D: Proto 


Gas: jWiD=4.9 


U: All 




56% 


23% 


44% 


2 % 


U: Starless 




68% 


40% 


89% 


54% 


U: Proto 




53% 


54% 


13% 


1 % 


Gas: A^iD = 


=4.9 


14% 


14 % 


56% 




Note. — 


D = 


driven, 


U = undriven 







obtain velocity dispersions and centroid velocities for 199 
prcstcUar and protostellar cores using NH3 (2,2), NH3 
(1,1) and C2S (2,1). They adopt a dust temperature of 
11 K, which is slightly lower than the assumed tempera- 
ture of 15 K used by K07. In comparison to p Opliiuchus, 
the Perseus star-forming region is much larger, 5 pc x 25 
pc, resembles a long chain of clumps with typical column 
densities of N-^^ ^ 3 x 10^^ cm~^, and contains a total 
mass of - 18,500 (Kirk et al. 2006). Using the total 
mass and assuming a cylindrical geometry {L = 25 pc 
and R = 2.5 pc) we obtain nn = 1-1 x 10^ cm~^ for 
Perseus, which we adopt for the simulation. We assume 
that Perseus is approximately in the plane of the sky; 
if it were randomly oriented then the expected value of 
the longest side of the cloud would be 50 pc. We take a 
temperature of 10 K for Perseus, since this is character- 
istic of the prestellar cores (R07). Equations (2) and (3) 
then imply that the simulation box has L = 2.8 pc and 
M = 825 M0, which is a relatively small piece of the 
total cloud. Since we are simulating only a small part 
of the Perseus cloud, we estimate the velocity dispersion 
in actual molecular gas from the average linewidth-size 
relation (eq. 6 for L = 5 pc), which gives tr = 1.1 km 
s~^ and Mid = 5.9. In comparison, our simulation box 
scaled to the Perseus average number density is less tur- 
bulent and only half the length of the shorter dimension. 
This difference in Mach number and cloud side yields a 
virial parameter for Perseus of a ~ 1, which is about 
60% of the value of our simulation box. 

5.2. Optical Depths 

In our analysis we make the assumption that the line 

transitions are optically thin. This approximation is ob- 
servationally validated for both the N2H+ and NH3 tran- 
sitions. For example, according to K07 the total optical 
depth, Ttot ~ 0.1 — 13, where rtot is the sum of the op- 
tical depths for each hyperfine transition. Thus, the av- 
erage optical depth for a given N2H+ hyperfine line is 
T = ''"tot/7 ^ 0.01 — 2, so that the majority of the lines 
are at least marginally optically thin. In particular, the 
isolated 101-012 hyperfine component used for velocity 
fitting has an optical depth of Ttot/9, and is therefore 
optically thin in all but the very densest cores. A07 re- 
port similar N2H+ total optical depths of Ttot ~ 0.1 — 30 
for p Ophiuchus. R07 find Ttot ^ 0.4 - 15 for NH3. The 
NH3 (1,1) complex has 18 hyperfine components so that 
most of the lines are at least marginally optically thin. 
For comparison, we report the total optical depth in our 
simulations for all three tracers in Table 4. We derive the 
optical depth for a given line by solving for the level pop- 
ulations as described in §3. Once these are known, the 
opacity in each cell for photons emitted in the transition 



(15) 

47r(wi - vo)vij ' 

where n, v, and X are the number density, velocity, and 
molecular abundance in the cell, Bji and Vij are the Ein- 
stein absorption coefficient and frequency of the transi- 
tion, and the observation is made in a channel centered 
at velocity Wobs that runs from velocity vi to Vq. The 
optical depth is given simply by computing this quantity 
in every cell, multiplying by the cell length to obtain the 
optical depth of that cell, and then summing over all cells 
along a given line of sight. As the table shows, for the 
most part the average hyperfine transition is optically 
thin in all tracers. The main exception is cores traced 
by NH3 in p Ophiuchus, which is marginally optically 
thick. As a result, we do not present results for NII3 
using the higher density p Ophiuchus scaling; the core 
velocity dispersion maps in Figures 5-8 are normalized 
to Perseus. 

In all other cases even the strongest hyperfine compo- 
nents have optical depths of order unity, and comparison 
with more detailed radiative transfer modeling than we 
perform indicates this is unlikely to significantly affect 
our results. Tafalla et al. (2002) model the emission and 
transfer of the same N2H"'" and NH3 lines that we use in a 
sample of starless cores in Taurus and Perseus whose con- 
ditions are similar to those produced by our simulations. 
They study the effect of the interplay between hyper- 
fine splitting and radiative trapping by analyzing the two 
limiting cases of negligible radiative trapping (whic:li we 
assume) and neglect of hyperfine splitting (which maxi- 
mizes radiative trapping). They find that the difference 
in the level populations they compute under these two 
assumptions is only a few tens of percent, a level of er- 
ror comparable to that introduced by uncertainties in 
the collision rate coefficients. We expect the errors in- 
troduced by our optically thin assumption to be compa- 
rable. 

5.3. Comparison of Second Moments 

Observationally, the second moments of cores are pre- 
dominantly subsonic in MCs, apparently independent of 
the amount of turbulence. For example, A07, measur- 
ing second moments in p Ophiuchus, find all values are 
smaller than 2cs with an average ctnt/cs = 0.5. Likewise, 
K07 report similar measurements for cores observed in 
Perseus, finding an average of <tnt/cs= 0.7 with a maxi- 
mum value of 1.7. Both our simulations find marginally 
sub-sonic distributions of second moments with slightly 
larger means than the observations (see Table 1). In com- 
parison, protostellar cores are observed to have a some- 
what broader distribution of second moments. K07 find 
that the protostellar cores in Perseus have a mean second 
moment of l.lcg and a maximum of 2.3cs. The protostel- 
lar objects that we observe in our driven simulation tend 
to have transonic second moments while in the decaying 
simulation they are supersonic. 

We use a KS test to compare the distribution of sec- 
ond moments for each of the simulation core populations 
with the observed core populations. We give the results 
in Table 5. Note that the A07 sample is comprised of 
only prestellar cores, while R07 observe both starless and 
protostellar cores but do not distinguish between them. 
Figure 11 shows the cumulative distribution functions of 
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TABLE 3 

Dispersion of first moments (centroid velocities) normalized to the large-scale gas dispersion. 
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Fig. 11. — Cumulative distribution function showing the total fraction / of cores with second moments, ctnt, less than or equal to the 
X coordinate value for simulated observations of p Ophiuchus and Perseus in N2H+ and NH3. The legends indicate by first letter whether 
the distribution is taken from K07, A07, R07, Undriven simulation, or Driven simulation. The tracer is also indicated when two different 
tracers are used. 



TABLE 4 

Total optical depth r through core centers for each 
normalization and simulated racer. 







Perseus 




p Ophiuchus* 




Ttot'' 


median 


min 


max 


median 


min 


max 


CI8O 


0.51 


0.08 


2.46 


0.35 


0.14 


1.05 


N2H + 


0.71 


0.07 


8.91 


7.27 


1.72 


29.44 


NH3 


8.37 


0.10 


63.49 


46.59 


10.61 


228.73 



^ Optical depths are reported for the distribution of starless cores 
only. 

''rtot is the sum of the optical depths through line center for each 
liypcrfine transition. For N2FI+ and NH3 with 7 and 18 hyperfine 
transitions, respectively, the optical depth is significantly reduced 
and generally optically thin for individual transitions. 

the core populations for some of the simulations and ob- 
servations. Although the medians of some of the second- 
moment distributions are fairly similar, KS tests of the 
core populations show significant disagreement in some 
cases. Overall, the distribution of second moments for 
the driven run is closer to observations of Perseus, while 
the decaying run is a better match for the p Opiuchus 
prestellar second moments. 
The physical origin of the poor agreement between the 



TABLE 5 

KS STATISTICS FOR THE DRIVEN AND DECAYING CORE SECOND 
MOMENTS (velocity DISPERSIONS) COMPARED TO THE 
OBSERVATIONAL COLLECTIONS OF CORES USING THE APPROPRIATE 
CLOUD NORMALIZATION AND SIMULATED TRACER. 
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simulations and observations appears to be that the sim- 
ulated protostellar second-moment distributions in either 
case do not have sufficiently narrow peaks. The pro- 
tostellar cores in the simulations are at the centers of 
regions of supersonic infall, which contradicts the obser- 
vations that show at most transonic contraction. Al- 
though the decaying simulation has a larger population 
of high dispersion protostellar cores, both simulations 
show almost equally bad agreement with the observa- 
tions. Tilley & Pudritz (2004), performing smaller decay- 
ing turbulent cloud simulations at lower resolution with 
self-gravity, analyze the linewidths of their cores using a 



15 



TABLE 6 

KS STATISTICS FOR THE DRIVEN AND DECAYING CORE FIRST 
MOMENTS (CENTROID VELOCITIES) COMPARED TO THE 
OBSERVATIONAL COLLECTIONS OF CORES USING THE APPROPRIATE 
CLOUD NORMALIZATION AND SIMULATED TRACER. 



Sample 


Cloud 


D 


U 


Starless 


p Ophiuchus (AO?) 


0.5 % 


6% 




Perseus {K07) 


48% 


12% 


Protostellar 


Perseus {K07) 


6% 


85% 


All 


Perseus (K07) 


0.8% 


7% 




Perseus (R07) 


7% 


3% 



similar simple chemical mode. They also find a number 
of cores with greater than sonic central linewidths. There 
are two possibilities for the discrepancy between the ob- 
served protostellar cores in our simulation and those ob- 
served in Perseus. In reality, forming stars are accompa- 
nied by strong outflows that may eject a large amount 
of mass from the core, leading to efficiency factors be- 
tween Ecore =0.25-0.75 (Maztncr & McKee 2000). Such 
outflows limit the mass of the forming protostar by this 
amount. Since we do not include outflows we naturally 
expect our sink particles to overestimate the forming pro- 
tostar mass by this factor and hence the maximum infall 
velocity, characterized by the second moment through 
core center. If we adopt a sink particle mass correction of 
3 (Alves et al. 2007), then the infall velocity will decrease 
by a factor of v^. This correction substantially reduces 
the number of protostellar cores with supersonic second 
moments from 53% and 70% to 23% and 39% for cores in 
the driven and decaying simulations, respectively. This 
correction brings the driven core sample closer in agree- 
ment with those measured by R07 and K07. A second 
possibility for the higher second moments is the lack of 
magnetic fields in oiir simulations. Magnetic pressure 
support could also retard collapse and decrease the mag- 
nitude of the infall velocities. However, the importance 
of magnetic effects is difficult to assess without further 
simulations. 

5.4. Comparison of First Moments 

In contrast, we find better agreement between simu- 
lations and observations for bulk core motions. When 
comparing the distributions of first moments, we first 
subtract out any large gradients in the sample as dis- 
cussed in §4.3. This is particularly important when com- 
paring to a large elongated cloud such as Perseus. We 
then shift the distributions so that median centroid ve- 
locity falls at and normalize the distribution to the bulk 
gas dispersion. For Perseus, we infer the bulk gas velocity 
dispersion for our simulation cr = 1.1 km s~^ by assuming 
the cloud falls on the linewidth-size relation and satisfies 
equation (6) with L equal to the transverse size of the 
cloud. For p Ophiuchus, we adopt the ^''CO line velocity 
dispersion of cr = 1.06 km s~^ (Loren 1989b). 

In Table 6, we report the KS agreement for the first 
moments of the observations and simulations. Since the 
simulations themselves are statistically similar to one an- 
other, both of the first moment distributions generally 
either agree or disagree with the observed population. 
Except in the case of the N2H+ driven data for p Ophi- 
uchus and the NII3 decaying data, the velocity-corrected 
data are fairly statistically similar to the observations. 
This suggests that the first-moment distributions do not 



strongly depend upon the details of the turbulence. In 
figure 12, we have plotted the cumulative distribution 
function of some of the first-moment distributions for 
comparison. The net core distributions show substantial 
overlap for both simulations and observational regions. 
The main source of disagreement with observations is the 
generally larger dispersions of the first moments in the 
simulations. In particular, the dispersion of the prestellar 
core first moments is a factor of ^ 2 larger than the that 
found by A07 in p Ophiuchus. However, because in both 
simulations the core-to-core velocity dispersion is smaller 
than the virial velocity of the cloud on large scales, we 
conclude that a sub- virial dispersion of first moments is 
not necessarily an indicator of global collapse. 

In some cases, the direct dispersion of the gas may be 
poorly observationally constrained and so a virial argu- 
ment is used to infer the gas dispersion. We find that 
normalizing the distributions to the virial gas disper- 
sion rather than the measured gas dispersion produces 
a significantly different result for the decaying simula- 
tion. Since the cloud gas is becoming more quiescent with 
time, the actual gas dispersion is sub-virial at late times. 
Thus, relative to the virial gas dispersion the decaying 
dispersion of first moments appears twice as sub-virial. 

6. DISCUSSION AND CONCLUSIONS 

We use isothermal AMR simulations to investigate the 
kinematics of cores in environments with and without 

driven turbulence. We simulate observations of these 
cores in the tracers C^^O, N2H+, and NH3 for the 
star-forming regions p Ophiuchus, 125 pc distant, and 
Perseus, 260 pc distant, with beam sizes of 26" and 31" , 
respectively. From the differences between cores in the 
two environments and in conjunction with observational 
results, we are able to draw a number of important con- 
clusions, some of which are relevant for observationally 
distinguishing between driven and decaying turbulence 
in star-forming clouds. 

We find that in both simulated environments the 
prestellar second-moment distribution is fairly narrow 
and peaked about the sound speed. Significant broad- 
ness of the protostellar second moment distributions is 
due to strong infall, such that many cores have central 
dispersions exceeding 2cs. Despite these commonalities, 
a KS test indicates that the driven and decaying prestel- 
lar and driven and decaying protostellar populations are 
dissimilar to one another. In contrast to the second mo- 
ments, a KS test indicates that the first-moment distri- 
butions in the two environments have some overlap: 13% 
confidence for protostellar cores and 44% confidence for 
prestellar cores. This similarity is an indication that the 
bulk core advection is decoupled from the gas motions 
inside the core. The similarity of the KS tests suggests 
that core first moments are not a good method for dis- 
tinguishing the two environments. 

Examining the gas dispersion in the core neighbor- 
hoods reveals interesting differences in the two simula- 
tions. We find that by the end of a global frcefall time the 
averaged velocity dispersion increases strongly towards 
the core center for decaying protostellar cores. How- 
ever, for decaying prestellar cores and all driven cores 
this trend is fairly flat or slightly increasing. Thus for 
both phases the driven cores are coherent, similar to ob- 
served cores (Kirk et al. 2007; Barranco & Goodman 
1998; Goodman et al. 1998), while the supersonic veloc- 
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Fig. 12. — Cumulative distribution function showing the total fraction / of cores with first moments, Vient, less than or equal to the 
X coordinate value for simulated observations of p Ophiuchus and Perseus in N2H+ and NH3. Each line is normalized to the appropriate 
large-scale gas dispersion, Vg, cither as measured (simulations) or as derived from the linewidth-size relation in equation (6). The legend 

format is similar to figure 11. 

ities observed in decaying protostellar cores are incon- 
sistent with observations. Thus, investigating the radial 
dispersion of protostellar cores may make it possible to 
discriminate between clouds with and without active tur- 
bulent energy injection. 

We find that the majority of the combined prestellar 
and protostellar distribution of second moments through 
the core centers for both environments are below 2cs, 
which agrees with the results of A07 and K07. How- 
ever, neither prestellar core distribution shows a signifi- 
cant confidence level of agreement with the observations. 

As shown in Table 5, we obtain sub-virial dispersions 
of the first moments for both total core populations like 
A07, however our core-to-corc dispersions arc approxi- 
mately a factor of 2 closer to virial. Although both runs 
produce sub-virial core-to-core dispersions, we have not 
shown that cither driven turbulence or the small virial 
parameter of decaying turbulence can produce avir as 
small as that found by A07. 

One interesting finding is that the protostellar cores in 
the decaying run have a core-to-core dispersion that is 
higher than the gas dispersion measured after a free-fall 
time. This is a result of the significantly larger dispersion 
of the protostellar cores compared to the prestellar cores, 
which may be a result of either increased scattering or 
of memory of the natal higher dispersion gas. This is in 
contrast to the driven prestellar cores, which have nearly 
the same dispersion as the gas, and the driven proto- 
stellar cores, which have a sub-virial dispersion. Thus, 
comparing the starless and protostellar core first-moment 
dispersion to the net gas dispersion is potentially a means 
for distinguishing the two environments. 

An effect that we cannot rule out is the importance 
of magnetic fields, which we do not treat in our simula- 
tions. In addition to seeding the initial clump mass spec- 
trum, the turbulence in our simulations provides sup- 
port against the cloud's self-gravity, a role that could be 



filled by either sustained turbulence or magnetic fields or 
both. The very small number of cores observed with 
supersonic second moments indicates that these cores 
are collapsing very slowly, a condition that we find is 
promoted by turbulent support but not throughout the 
entire core collapse process. At present, little compu- 
tational work has been done to study line profiles for 
turbulent cores with magnetic fields. Tilley & Pudritz 
(2007) present central line profiles for a few cores formed 
in self-gravitating magneto-hydrodynamic cloud simula- 
tions but do not have many statistics. Our simulations 
also neglect protostellar outfiows, which may have an 
effect on the total core mass and hence the velocity dis- 
persion of the infalling gas in the core center. 

Another possible source of the quantitative disagree- 
ment between observations and our simulations is geom- 
etry. Periodic boundary conditions may do a poor job 
representing whole, pressure confined molecular clouds. 
Certainly, the star-forming region of Perseus is more fil- 
amentary than round. Further, the cloud Mach numbers 
for both regions are somewhat uncertain, and it may be 
necessary to match the Mach number of the simulation to 
the cloud more exactly to get better quantitative agree- 
ment. 

Overall, we find that the driven simulation agrees bet- 
ter with the cores in Perseus, while the decaying simu- 
lation agrees slightly better with the pre-stellar cores in 
p Ophiuchus (our data do not include protostellar cores 
there). Our results indicate that the decaying simida- 
tion produces a population of protostellar cores with su- 
personic velocity dispersions that is largely inconsistent 
with the observations of protostellar cores in Perseus. To 
reach a firmer conclusion on the validity of driven or de- 
caying turbulence will require more complete data on a 
larger sample of clouds as well as simulations that allow 
for magnetic fields, outflows, and thermal feedback from 
the protostars. 
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APPENDIX 



STATISTICAL EQUILIBRIUM FOR MOLECULES WITH HYPERFINE STRUCTURE 

As discussed in Tafalla et al. (2002) and Keto et al. (2004), hyperfine splitting in a molecule introduces two compli- 
cations on top of the standard calculation of statistical equilibrium. First, hyperfine splitting of a transition reduces 
its optical depth by breaking the line into multiple components. The frequency separation between the components 
means that photons generated by a transition from level ia to level j(3, where the Roman index refers to the parent 
level and the Greek to its hyperfine sublevel, have a reduced probability of being resonantly absorbed by molecules in 
state i that are not in hyperfine sublevel /3. Under our assumption that all compcments mc optically thin, however, 
this effect is not significant. We discuss the extent to which this approximation holds, and how our results might be 
modified in cases where it fails, in § 5.2. 
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A second, practical complication is that collision rate coefficients between different hyperfine sublevels are generally 
unknown. Only the total rate coefficients summing over all hyperfine states are known. This makes it impossible 
to perform a true statistical equilibrium calculation without introducing additional assumptions, the most common 
of which is that the individual hyperfine sublevels are simply populated in proportion to their statistical weights. 
Observations along some sightlines show that this approximation generally holds for NH3 and that deviations from it 
for N2H+ are only of order 10% (Tafalla et al. 2002; Keto et al. 2004). 

Under the assumption of an optically thin gas, the equation of statistical equilibrium for a molecular species with 
hyperfine structure is 

^("■Ha^j/Sia + ^j/3ia + BjjSialcMB) fj^ 
3 



k (3 



1, 



(Al) 
(A2) 



where a set of four subscripts iajfi indicates a transition from state i, hyperfine sublevel a to state j, hyperfine sublevel 
j3. The assumption that the hyperfine sublevels are populated in proportion to their statistical weight then enables us 
to write ^ 

fia — fit (■^"^) 
9i 

where gi^ is the statistical weight of sublevel ia, gi — gia is the total statistical weight of all hyperfine sublevels 
of level i, and fi = is the fraction of molecules in any of the hyperfine sublevels of level i. If we make this 

substitution in equations (Al) and (A2), then they become 
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If the hyperfine sublevels of state i are populated in proportion to their statistical weight, then the total transition 
rate from all hyperfine sublevels of state i to any of the sublevels of state j are given by 
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Now note that (A4) represents one independent equation for each state i and each of its hyperfine sublevels a. If we 
fix i and add the equations for each hyperfine sublevel a, then equation (A4) simply reduces to 



X^(?iH2ftfc + ^ik + Biklcvm) 



fi 



(A9) 



the same as the equation for an optically thin molecule without hyperfine splitting, provided that the rate coefficients 
are understood to be summed over all hyperfine sublevels. 



